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ABSTRACT 

Photo-electric absorption of blackbody photons is an important process which hm- 
its the acceleration of ions under the space-charge limited flow boundary condition 
at the polar caps of pulsars with positive corotational charge density. Photo-electric 
cross-sections in high magnetic fields have been found for the geometrical conditions 
of the problem, and ion transition rates calculated as functions of the surface temper- 
atures on both the polar cap and the general neutron-star surface. The general surface 
temperature is the more important and, unless it is below 10^ K, limits the accelera- 
tion electric field in the open magnetosphere to values far below those needed either 
for electron-positron pair creation or slot-gap X-ray sources. But such ion beams are 
unstable against growth of a quasi-longitudinal Langmuir mode at rates that can be 
observationally significant as a source of coherent radio emission. 
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1 INTRODUCTION 

Two previous papers have given the results of a study of 
reverse-electron flow in isolated neutron stars with posi- 
tive polar-cap corotational charge density pa J, in which 
free emission of ions maintains the boundary condition for 
space-charge limited flow. The sign of this quantity depends 
on the relation between rotation spin il and polar mag- 
netic flux density B. We refer to Goldreich & Julian (1969) 
for its Euclidean-space definition, and to Muslimov & Tsy- 
gan (1992) and Muslimov & Harding (1997) for the general- 
relativistic modification used here. Most publications have 
assumed the parallel case, • B > 0, with pcj < and 
electron acceleration. But there does not appear to be any 
reason why the antiparallel case with pcj > and ion and 
positron acceleration should not be present in the neutron- 
star population. The two earlier papers (Jones 2010a, 2011; 
hereafter Papers I and II) considered in some detail the con- 
sequences of proton formation in the electron-positron show- 
ers produced by the reverse-electron flux at the polar cap 
and its relation with observed radio-pulsar phenomena. But 
an important omission was the effect of reverse-electron flow 
on the electrostatics of particle acceleration under the space- 
charge limited boundary condition, and this is the subject 
of the present paper. 

In this paper, the subscripts || and _L are with refer- 
ence to the local magnetic fleld direction. Under the surface 
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electric-fleld boundary condition 7^ assumed in the 
classic paper of Ruderman & Sutherland (1975), reverse- 
electron flow does not modify ion or positron acceleration in 
any important way because the surface value of the accel- 
eration field E|| is unconstrained. However, the presence of 
a reverse-electron negative charge density under the surface 
boundary condition Ey = of space-charge limited flow can, 
in certain cases, reduce the acceleration field to values so low 
that there is no possibility of electron-positron pair creation 
within the open magnetosphere above the polar cap or of 
a slot-gap X-ray source. Photo-electric absorption of polar- 
cap and neutron-star surface blackbody photons by acceler- 
ated ions is the source of the reverse-electron flux and the 
process is analogous with electron-positron pair creation in 
its effect on the acceleration fleld. Heating of the polar-cap 
surface by a reverse flux of photo-electrons has been con- 
sidered previously by Cheng & Ruderman (1977), also by 
Muslimov & Tsygan (1987) in relation to an instability of 
a liquid surface under the boundary condition Ey 7^ 0. The 
present paper is concerned not with heating effects but with 
the electrostatic consequences of the negative charge density 
that reverse-electron flow implies. 

There have been many calculations of the photo-electric 
cross-section for atomic hydrogen in the high magnetic fields 
and temperatures of model neutron-star atmospheres. We 
refer to Potekhin & Pavlov (1997) for a critical survey of 
this work and to Potekhin, Pavlov & Ventura (1997) for 
a study of the adiabatic approximation used by many au- 
thors. Although model atmospheres containing low-Z ele- 



2 P. B. Jones 



ments have been constructed (see, for example, Mori & Ho 
2007), there appear to be no published explicit photo-electric 
cross-sections for such elements in high magnetic fields. Also, 
the intervals of electron separation and photon energies rel- 
evant to the accelerated ions in the present work differ from 
those of neutron-star atmospheres. But the requirements of 
this paper have some simple features compared with the lat- 
ter problem. It is not necessary to consider the complexities 
of thermal ion motion transverse to the field and, in par- 
ticular, the photon momentum component k± is such that 
in relation to the electron cyclotron radius rs, the condi- 
tion k±rB <C 1 is always satisfied. Hence simple expressions, 
valid only in this limit, have been derived and arc given in 
Appendix A. The zero-field cross-sections have been com- 
puted extensively (see, for example, the work of Reilman 
& Manson 1979) and in approximations much superior to 
those of Appendix A, but nevertheless we show that use of 
the high-magnetic field expressions is essential in this paper. 



The region of open magnetosphere above the polar cap 
is assumed to have a circular cross-section of radius uo{z), 
where z + _R is the coordinate on the magnetic axis and R is 
the neutron-star radius. The acceleration field has two com- 
ponents. The first, present at very low altitudes z <^ wo(0) 
above the polar-cap surface, is the inertially-generated field 
described by Michel (1974). The second, and much the more 
important, contribution is produced by the Lense-Thirring 
effect modification to the corotational charge density first 
recognized by Muslimov & Tsygan (1992) and extends over 
moderately large distances z ^ R. The computational re- 
sults of this paper are contained in Table 1, which gives the 
charge evolution of typical ions as they are accelerated. 



The photon source is modelled as a polar cap of ra- 
dius Mo(0) and proper-frame temperature Tpc embedded in 
a spherical neutron-star surface of uniform temperature Ts. 
We find that the latter temperature is much the more sig- 
nificant but, unfortunately, it is essentially unknown in the 
case of the general isolated neutron-star population. Black- 
body emission has been detected reliably in only a small 
number of pulsars with ages less than 10® yr (see the review 
of Yakovlev & Pethick 2004) and even then there remains 
the problem of deciding whether the source is the polar cap 
or the whole surface. 



Section 2 considers the simplest possible acceleration 
process: one that is time- independent and in a common state 
of plasma composition over the whole polar cap. It is shown 
that, in this case, acceleration is inadequate to produce sec- 
ondary electron-positron cascades except for values of Tg 
that are too small to be observable. The relation between 
this result and earlier work on space-charge limited flow in 
the f2 ■ B < case (Papers I and II) is discussed in sec- 
tion 3. The flux of ions that are relativistic but not ultra- 
relativistic is unstable against growth of a quasi-longitudinal 
Langmuir mode at rates that appear large enough to pro- 
duce Langmuir solitons and strong turbulence (Weatherall 
1997, 1998; Melrose 2000; Melikidze Gil & Pataraya 2000; 
Asseo & Porzio 2006). 



2 ACCELERATION AND PHOTO-ELECTRIC 
ABSORPTION 

Solutions for the charge density p, for pcj and the elec- 
trostatic potential 3> have been obtained for space-charge 
limited flow in a dipole field by Muslimov & Harding (1997) 
and by Harding & Muslimov (2001). Those authors consid- 
ered electrons, for which the inertial component of the ac- 
celeration field is negligible. For ion acceleration, the inertial 
component, though less important than the Lense-Thirring 
induced term should be included. The expressions given by 
Muslimov & Harding arc complicated and therefore wo shall 
make use of the boundary condition and geometrical form 
of the open fiux-line magnetosphere to obtain a very sim- 
ple approximate expression for # which is adequate for the 
purposes of this paper. We calculate the photo-electric tran- 
sition rate for an ion accelerated to Lorentz factor 7 at alti- 
tude z above the magnetic pole. 

There is little information about electron separation 
energies for multiply ionized atoms in the magnetic fields, 
lO^'^ < B < 10^'^ G of interest here. The most recent calcu- 
lations (Medin & Lai 2006) are restricted to atomic rmmbors 
Z = 6 and 26 and further, in the latter case to i? ^ 10^* G. 
We therefore adopt the expression used previously (Jones 
1981) for the separation energy ei, 



ei = 0.35B 



1/2 



Z-l 



keV, 



(1) 



in which I = 0,1,2,.. is the electron spatial degeneracy 
quantum number. When summed over all /, equation (1) 
gives agreement with the directly computed total energies 
of atoms with atomic numbers Z = 10 and 20 (see Jones 
1985) that is adequate for the work of this paper. 

The inertial component of the acceleration field can be 
found at small z <^ wo(0) by means of a one-dimensional 
approximation to the Poisson equation which neglects the 
boundary condition on the cylindrical surface u = uo separ 
rating open and closed flux lines. 



Ze 



1/2 



(2) 



for ions of mass M and charge Z. The field increases from 
the boundary condition Ey = to this constant value in a 
very small distance 



Zi 



\\%-KpGjZe) 



(3) 



Photo-electric absorption of blackbody photons modifies 
by introducing the negative charge density of the electrons 
accelerated toward the polar-cap surface. To estimate their 
effect, we can assume that the transitions all occur at a 
unique altitude z = h, increase the ion charge from Z to Zh, 
and reduce the field to -Bii = at « = h. Then, 



and the ion Lorentz factor at a = /i is, 
7ft « 1 + o 



(4) 



(5) 



'^Zh-Z 

The electron chemical potential in the upper, low-density. 
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regions of the atmosphere at the polar cap surface is, 
He = fesT In nrgNe 



rnk B T 



(6) 



at electron number density A^e. We have estimated Z by 
comparing the value of fi^ on the surface of last scattering 
in the atmosphere prior to acceleration with the e; given by 
equation (1). On this basis, and because the dependence of 
He on B is only logarithmic, we have adopted the condition 
£( = 26kBT to obtain the Z used in Table 1. The polar-cap 
radius is that given by Harding & Muslimov (2001), 

M0)-(-SM^X'\ (7) 



cPf(l))J ' 

in which P is the rotation period and /(I) = 1.368 for a 
neutron star of mass IAMq and radius i? = 1.2 x 10® cm. 

The inertial component of E\i ceases to be significant at 
h ~ 10'' cm. In principle, calculation of the electric potential 
$ at z > /i requires the solutions given by Harding and 
Muslimov. But at z ^ mo(0) a very simple expression can 
be found on the basis that the open flux lines are contained 
within a narrow cylindrical tube whose cross-sectional area 
increases only slowly with z. Given the boundary condition 
$ = on the surface u = uo and that the radial electric field 
is E± ^ £|| , application of Gauss's theorem to a section of 
the tube at altitude z gives the approximation, 

$ = TT (wg (Z) - U^) (piz) - PC J (Z)) , (8) 

in terms of the difference between the charge density p{z) 
and the Goldreich- Julian charge density pgj{z). We refer 
to Harding & Muslimov (2001) for this quantity at polar 
coordinate r = r]R and adopt the expression. 



p - pa J = 



cPar]^f{l) 



COS'!/', 



(9) 



in which we have neglected a second term containing sin-ip 
which is small at altitudes r] of order unity. Equation (9) 
assumes a dipole field with surface value B. Here, k is the 
dimensionless Lense-Thirring factor and tp is the angle be- 
tween $1 and B. This expression is a satisfactory approxi- 
mation at « > ito(O) and z < 3R. At these altitudes, the 
function /(j?)/a/(l) including the red-shift factor a can be 
well-approximated by unity and is a slowly varying function 
of rj. 

Equations (8) and (9) define the acceleration of an ion 
and give the relation between the altitude rj and the Lorentz 
factor 7 and we restrict transition rate calculations to points 
on the magnetic axis (w = 0). The neutron-star surface is 
assumed to be Lambertian. Then unit element of surface 
area at a point which subtends an angle with the magnetic 
polar axis produces a photon number flux. 



Jo{oJo,6) = 



uJon{ujo) cos(C + 0) 



(10) 



47r3c2J?2 1 -h 772 -277 cos C' 

at altitude z — r — R on the axis, in which the photon 
occupation number is 'n{ujci) for photons of angular frequency 
iJo- The photon momenta from this element are at an angle 6 
with the magnetic a:xis. Lorentz transformation parallel with 
B to the rest frame of an ion on the axis moving outward 
with Lorentz factor 7 gives the flux in that frame, 



where oj and x replace uio amd 9, and /3 is the ion velocity 
in units of c. The photo-electric transition rate transformed 
back to the neutron-star frame at an altitude r? — 1 on the 
magnetic axis is. 



7/ 



27ril'sinCdC / J(a;, x)a(B, w, x), (12) 
Jo 

in which the upper limit of integration is given by cos (^max = 
1/r] and the integral is in two sectors; the polar cap with tem- 
perature Tpc and radius given by equation (7), and the re- 
maining visible neutron-star surface of temperature Ts . The 
cross-section is cr = (ag + (Ji + o-\)/2. Here, the subscripts 
denote the Landau quantum number of the final electron. 
The photon polarization is denoted by the superscripts o, h, 
respectively perpendicular to or parallel with kx B. General- 
relativistic corrections to photon propagation, including the 
red-shift factor, have been neglected here in equations (10) 
- (12). Equation (12) gives the transition rate as a function 
of the coordinate r? on the magnetic axis and the Lorentz 
factor 7 of the ion. Owing to the geometrical restriction of 
equation (12) to points on the magnetic axis, transition rates 
for ions accelerated off-axis at finite u{z) would be in error 
to a modest extent for the polar-cap contribution, but for 
the much more significant whole-surface term, the error is 
negligible. We refer to Appendix A for details of the indi- 
vidual cross-sections. Approximations and sources of error 
are also discussed there, as is the comparison with zero-field 
cross-sections. 

The integration in equation (12) is over the polar cap 
with temperature Tpc and the remaining visible part of the 
neutron-star surface at temperature Ts- From the transition 
rates and the relations between 7 and 2 given by equations 
(2) - (5), (8) and (9), two optical depths have been deter- 
mined for three values each of nuclear charge and polar-cap 
magnetic flux density. These are. 



Vi{z)dz 



(13) 



which defines Z^ for a given Z. Its dependence on Ts is 
negligibly small principally because the factor cos(C-|-^) ^ 0. 
The second optical depth is. 



-/" 



Vi{z)dz 



(14) 



which gives Zoo and is almost completely independent of 
Tpc owing to the unfavourable nature of the Lorentz trans- 
formation for polar-cap photons &t z^ uo (0) . In both cases, 
integer charges are determined by whether or not the optical 
depth concerned is greater or less than unity. This is usu- 
ally obvious because optical depths decrease rapidly with 
increasing ion charge and electron separation energy. 

Equations (4) and (5) are based on the effect of reverse- 
electron negative charge density at 2: <IC mo(0). Photo-electric 
absorption at « ^ wo(0) also introduces a reverse-electron 
flux and negative charge density in any interval of z in which 
there is a finite E\\. Thus puo for ultrarrelativistic ions is 
not constant but is a function of z determined by the local 
transition rate, 



J(w,x) =7Jo(wo,6>)(l-/3cos6>), 



(11) 



1 d{pul) _ SF; 
pul dz ^ cZ{z) 
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and so reduces the acceleration field obtained from equation 
(8). There can be acceleration in any interval of z only if 
photo-electric rates are such that d$/dz < and this con- 
dition leads at once to the inequality, 



E 



2r, 

cZ{z) 



dz 



d ( ^ 1 
</t— 1-^ 



77 



(15) 



in which the fractional error in the right-hand term is of 
order k relative to unity, and Z{z) is the local value of the 
ion charge. Because the principal aim of this calculation is to 
see if there can be sufficient acceleration to produce electron- 
positron pair creation, wc adopt the following procedure. We 
scale $ so as to give an ion Lorentz factor 7 = 500 at Zmax, 
that is, a potential diflference of about 10* GeV. Then the 
lowest ion charge that has less than unit optical depth in the 
interval h < z < Zmax is defined as Zoo. But if this acceler- 
ation is to be possible, photo-electric rates must not exceed 
the limit defined above. On integrating (15), we find that Zoo 
must satisfy the approximate inequality Zoo — Zh < kZoo/2 
for consistency. This condition is maintained by an equi- 
librium such that the acceleration field satisfying Poisson's 
equation adjusts to those values that give the ion acceler- 
ation and photo-electric transition rates necessary for self- 
consistency. To enlarge on this, we can note that if photo- 
ionization were complete at an altitude Zc, acceleration of 
the bare nucleus would continue a.t z > Zc, the potential dif- 
ference remaining for this being approximately ^{z) — ^{zc)- 
But integration of (15) for this case gives the approximate 
constraint. 



so that values of Z — Zh for which continued acceleration is 
possible are severely constrained. The equilibrium wo have 
mentioned above means that the acceleration field given by 
Poisson's equation adjusts at all z to that required for photo- 
ionization transition rates no larger than those satisfying 
the inequality (15). Reference to Table 1 shows that con- 
tinued acceleration over a potential difference approaching 
the maximum given by equations (8) and (9) is not possible 
except for low B and high values of Tpc. 

The optical depth integration is continued to Zmax = 3J? 
although the major contribution is from much lower alti- 
tudes. But the local magnetic flux density is much reduced 
near Zmax and therefore the separation energy used is the 
greater of that given by equation (1) or a zero-field expres- 
sion given by a Rydberg formula in terms of the screened 
nuclear charge and the appropriate principal quantum num- 
ber for the electron shell. 

The Table shows that Z^ differs little from Z for any 
values of B and Tpc. There are two reasons for this. Firstly, 
the cross-section for transitions to the lowest Landau state, 
the only final electron state accessible at low values of 7, con- 
tains the factor sin^ Xi which is always small. Thus use of the 
zero-field cross-section in this region would be quite incor- 
rect. Secondly, as previously explained in relation to equa- 
tions (4) and (5), the reverse-electron negative charge den- 
sity itself limits acceleration at z < h. Most photo-electric 
absorption occurs at z ^ «o(0), is Ts-dependent, and con- 
sists of transitions to the n = 1 Landau state whose cross- 



sections contain cos^ x- These, of course, require acceleration 
to higher Lorentz factors. 

The estimate of the Lense-Thirring parameter used by 
Muslimov & Harding (1997) is ft = 0.15. Values of Zoo are 
given for neutron-star surface temperature Ts = 1.0 x 10'"' 
K in column 6 and by reference to the integrated form of 
(15) show that Zoo — Z^ is such that acceleration is not al- 
ways possible even at this low temperature. At 2.0 x 10^ K, 
acceleration occurs only in the limited number of cases in 
which high values of Tpc give Z^ = Z. Because Z^ differs 
little from Z, we can see that the basic uncertainty is in the 
latter quantity. Estimates more refined than our condition 
derived directly from equations (1) and (6) might give lower 
or higher values of Z. In the former case, values of Zh would 
also be smaller, hence the possibility of significant acceler- 
ation would be reduced. In the latter. Zoo — Zh would be 
reduced. However, if photo-electric ionization is complete at 
z = h, the energy fiux of reverse electrons is very limited. 
Typically, for Zh — Z = 1, equation (5) shows that the re- 
verse electron energy input to the polar cap is no more than 
~ 20 GeV per ion, giving only low rates of pair creation by 
inverse Compton scattering or by (n,7) reactions. 

The values of polar-cap magnetic flux density consid- 
ered in Table 1 were chosen so as include the median and 
mean values of logj^g B for the Wang, Manchester & John- 
ston (2007) catalogue of pulsars that exhibit nulls. For fields 
much smaller than 10^^ G, the approximation k±rB -C 1 
made in Appendix A fails. Landau level spacings become 
small, and the photo-electric cross-section approaches the 
zero-field value. Thus transition rates remain high. On the 
basis of Table 1, we find, with some confidence, that the 
Muslimov- Tsygan acceleration field is in general so reduced 
by electron production that electron-positron pair creation 
is not possible. In this state, which we assume to be time- 
independent, there is an equilibrium between the photo- 
electric transition rates and the ion acceleration needed to 
produce them. In many cases, ion Lorentz factors at Zmax 
may differ little from 7/1. Our finding is limited, of course, to 
neutron stars whose whole-surface temperature is Ts ^ 10^ 
K. But a temperature of 10^ K is problematical from all 
points of view. It is unobservable in a compact source given 
present techniques and the energy flux needed to maintain it 
is small and easily supplied by many possible dissipative pro- 
cesses. Also, measured photon spectra do not always allow 
the unambiguous separation of polar cap and whole-surface 
radiation. 



3 CONCLUSIONS 

The results described in the previous Section are a further 
example of the diverse forms of polar-cap acceleration un- 
der the space-charge limited flow boundary condition that 
exist for the f2 ■ B < case. It was assumed that the com- 
position of the accelerated plasma is uniform over the whole 
polar-cap area. Before examining possible observable conse- 
quences of the uniform state, we need to mention briefly the 
impact of these results on the polar-cap acceleration which 
was considered at length in Paper II. 

We first comment on the assumed stability of the state 
described in Section 2. In the physical processes concerned, 
there is no obvious characteristic time apart from ion flight 
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Table 1. For a range of values of polar-cap magnetic flux density, 
temperature and nuclear charge Z, values of the ion charge in 
local thermodynamic equilibrium immediately before acceleration 
(Z) and at the end of the first stage of acceleration (Zf^) are 
given. The ion charge at the end of acceleration Z 

oo IS given in 

the right-hand four columns for arbitrarily selected values of the 
general surface temperature Ts subject to the explanation and 
qualifications discussed in Section 2. For Tg ^ 5.0 x 10^ K, the 
ion charge is always Zoo = Z. 
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times which are of the order of Rc~^ ~ 4 x 10~® s. This 
is many orders of magnitude shorter than the proton diffu- 
sion time constant Tp defined in Paper II. Thus we do not 
believe that the possibility of instability arising from photo- 
electric absorption can materially change the conclusions of 
that work. 

Decay of the giant-dipole resonance formed in reverse- 
electron showers produces protons, and it was shown in Pa- 
pers I and II that this process leads to temporal instabilities 
in the composition of the plasma accelerated. It was pro- 
posed that the consequences of these include the phenomena 
of nulls and sub-pulse drift observed in some radio pulsars. 
The number of protons created per unit ion charge accel- 
erated, defined as K in Papers I and II, is the significant 
parameter. It was shown in II that, in general, the state of 
the polar-cap must be time-dependent, with any given el- 
ement of area alternating between the acceleration of ions 
and positrons or of protons. In such an element, proton ac- 
celeration produces no reverse-electron flux and so must be 
of limited duration. It was assumed in Paper II that accel- 
eration would be through the full potential difference given 
here by equations (8) and (9), but this obviously needs some 
qualification. At any instant, there is positron and ion accel- 



eration from only a fraction (K + l)^^ of the whole polar-cap 
area. Thus the photo-electric production of negative charge 
is confined to this fraction and for K ~ 5—10 the charge den- 
sity averaged over the cross-sectional area -ivUq{z) does not 
differ much from that given by equation (9). Very roughly, 
the modified integrated form of (15), 



(16) 



is now the necessary condition for the existence of a non- 
zero acceleration potential difference. The actual value of 
Zoo — Zh determines the amount by which the equilibrium 
acceleration potential is reduced from the full potential $ . 
Inspection of Zoo — Zh values given in Table 1 shows that, 
for if ^ 5 and with the exception of Tpc = 1.0 x 10® K, the 
acceleration potential is usually a large fraction of $. 

The uniform state considered here in Section 2 is of 
interest in connection with the usual assumption that sec- 
ondary low-energy electron-positron pair creation is neces- 
sary for coherent radio omission. It has no obvious mecha- 
nism for electron-positron pair formation and the assump- 
tion might well be that, in an isolated neutron star, it would 
correspond with a long null state of radio and polar-cap X- 
ray emission. 

However, the possibility that cold ions might be im- 
portant in the growth of unstable quasi-longitudinal plasma 
modes was considered many years ago by Cheng & Ruder- 
mam (1980), but it appeared that the longitudinal effective 
meiss M'f^ of the accelerated ions was too large to give the 
necessary growth rates. However, the small Lorentz factors 
found here prompt reconsideration because the longitudi- 
nal effective mass is much reduced although the energy flux 
of the ions remains high in comparison with typical pulsar 
radio luminosities. 

In the uniform composition state, acceleration fields are 
everywhere small so that the reverse-electron energy fiux 
gives values K < 1 and a mixed plasma of ions and pro- 
tons for which the instability arguments of Papers I and 
II are not necessarily valid. The plasma then has two cold 
components, ions and protons with small Lorentz factors -,i 
and 72, respectively, with 71 « O.572, and obviously satisfies 
the relativistic Penrose condition (see Buschauer & Benford 
1977). But because the plasma is cold, we can easily find the 
growth rate of the quasi-longitudinal Langmuir mode stud- 
ied by Asseo, Pelletier & Sol (1990) for which the dispersion 
relation is, 



LO2 



(w-fc||/3l)2 {ui-kniS'zy 



1 2 

= k±, 



(17) 



for angular frequency oj, wave vector k, and at ion and pro- 
ton velocities fii and {32, respectively. The remaining quan- 
titles are w. 



7j '^^^LOi in which. 



2 AirNiqt 



Mi 



(18) 



where Ni, qi and Mi are respectively the number density, 

charge and mass of each component i = 1,2. The rest- 

1/2 

frame plasma frequency for a given component is 7^ coi 
and therefore, following Asseo et al, we consider the case 
k± = and = 2a;i7i, defining a new variable s through 
the relation. 



U! - k»l3i = + s) 



(19) 
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so that the dispersion relation becomes, 
1 C 



1 



(S + 1)2 (S + 







(20) 



in which fi — 71/72 and C = lo*^ jixi^ . Inspection of equa- 
tion (20) shows that this quartic has two real roots which 
can easily be found numerically thus giving the remaining 
complex roots. In general, the growth rate found from the 
complex roots is of the order of lO"'^^*, but is smaller if 
another ion species with qi/Mi 7^ q2/M2VS, substituted in 
place of the proton beam. Thus for a typical ion charge of 
q\ /e = 20 and a proton beam of number density N2 = 5Ni 
wc can choose C — 0.045 and /u = 0.2, for which the ampli- 
tude growth rate is, 



Imw = ±0.18wi = ±4.0 X 10^ 



1/2 



(21) 



rad s~^, at radius r, where S12 is the polar-cap magnetic flux 
density in units of 10^^ G and P is the rotation period. The 
amplitude growth factor, assuming constant 71 the interval 
< < Zmax, is expA, where. 



T""""' , Imw 
/ dz 

Jo ^ 



= 3.2 X lO'' 



B12 

rlP 



1/2 



\R+ Zmax ) ) 



(22) 



This is large for 71 < 10. Values of the mode frequency are 
relatively small, w w 271 w*, and may lie well below 400 
MHz, but it would be interesting to see if evidence for it is 
present in pulsar spectra at low frequencies. 

We have seen that the presence of ion and proton beams 
that are relativistic but not ultra-relativistic, previously un- 
expected in pulsars, can be associated with the growth of 
unstable quasi-longitudinal modes having growth rates high 
enough for the formation of Langmuir solitons and strong 
turbulence. We can refer to the review of Melrose (2000) 
and to the papers of Weatherall (1997) and (1998), Me- 
likidze et al (2000) and Asseo & Porzio (2006) to note that 
the collapse of such solitons is widely thought to be a plau- 
sible source for the coherent radio emission. It is usually 
assumed that the instability occurs in secondary electron 
and positron beams with Lorentz factors of the order of 10^ 
but it is the macroscopic fields of the solitons that lead to 
coherent emission independently of the particle type from 
which they are formed. 
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APPENDIX A: PHOTO-ELECTRIC 
ABSORPTION IN HIGH MAGNETIC FIELDS 

The presence of a high magnetic field greatly modifies the 
final states that are accessible to the electron in photo- 
electric absorption. In particular, the Landau level spacing, 
fibJB = 11.6 keV at 10^^ G, is large compared with black- 
body photon energies. In the rest-frame of the accelerated 
ion, the photon momentum k is almost parallel with B and 
this much reduces transition rates to the n = Landau state. 
We adopt the Johnson & Lippmann (1949) spinor solutions 
for a free electron in a uniform magnetic field. These are 
satisfactory here because we require only the cross-section 
(To to the unique n = ground state, and the sum of the 
cross-sections to the two degenerate n = 1 states. Photon 
polarization states are defined as those perpendicular to, or 
parallel with, k x B and arc denoted by superscripts (a, V) 
respectively. For the explicit form of the solutions and of the 
interaction with the electromagnetic field it will be conve- 
nient here to refer to Appendix A of Jones (2010b). 

Although with later comment, we assume the adia- 
batic approximation for the partially-ionized ground state 
in which the electrons occupy single-particle states each 
with Landau quantum number n = 0. Recoil is neglected 
so that in the ion rest frame, the photo-electric thresholds 
are e; -|- nfiUB , where ei is the separation energy given by 
equation (1) and n is here the electron final state Landau 
quantum number. Following equation (A6) of Jones (2010b), 
the transition matrix element from a photon with angular 
frequency w and polarization e"'** to the n = ground state 
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with electron longitudinal momentum p is, 
r±dr±dcj)e^-'-'''-'- 



*t_i,o(p)a • e"'''*_i,o(p - fc||)C?(p - fell) 



(Al) 



in cylindrical polar coordinates, where G is the Fourier 
transform of the model bound electron state v^exp(— 
in which ti^K^ = 2mei. The dependence of the cross-section 
on the charge of the ion is present solely through the Fourier 
transform. The model function, which replicates the correct 
asymptotic form of the true function, given in this Appendix 
by the parameter k, is an adequate approximation for the 
purposes of this paper. The Johnson & Lippmann spinor for 
Landau number n > is, 



*-l.n(p) = 







where p\ = 2mn?iu!B, and the normalized functions ipni are 
the non-relativistic radial solutions of the Schrodinger equa- 
tion for a free electron in a uniform magnetic field, in which 
I is the spatial degeneracy quantum number (sec equations 
A3 - A5 of Jones 2010b). The azimuthal angle-dependence is 
exp(— i(/ — n)(f)). The ion structure consists of spatially well- 
defined coaxial shells with sequential values I = 0,1, 2.... and 
energies given by equation (1). The possibility of a hole state 
is not considered because its lifetime would be many orders 
of magnitude smaller than the photo-electric transition rates 
calculated here. 

For the specific geometry of the photo-electric process 
considered here and for high fields B ^ 10^^ G, the exponent 
in the matrix element, which is of the order of fexrs is always 
very small compared with unity and is set equal to zero. This 
introduces the selection rule 5{n — l)=0 which leads to an 
easy evaluation of the transition rates. Hence, immediately 
from equation (Al), we obtain the cross-section. 



0-0 



27re^ sin^ X 
cpmuj 



K%2p- 



+ ■ 



K^(2p-|-fe|l 



(«;2-h(p-fe||)2)2 (K2-h(p-hfe||)2) 



(A2) 



for transitions to the n — Landau state for polarization 
perpendicular to k x B. In the parallel polarization case, the 
cross-section to the n = state is ctq = 0. Proceeding in a 
similar way, the cross-sections to the two degenerate n = 1 
states are found to be. 



2776^ cos^ X 

cpmuj 

K^ipl+kf.) 



+ 



K'{pi+kh 



(K2 + (p_fc||)2)2 (^2 + (p + fc||)2) 



(A3) 



and (Ti = CTl. It is to be emphasized again that these ex- 
pressions are valid only under the condition that fexrs <C 1 
which is valid for the photo-electric process considered here. 
These cross-sections, not too far above threshold and at 10^^ 
G, do not difi'er greatly from those given by the zero-field 
Kramers formula apart from the presence of the sin"^ x fac- 
tor in cTo, and the zero parallel-polarization cross-section ctq. 
These factors, within the adiabatic approximation, are a di- 



rect consequence of the relative orientation of the photon 
polarization vector and the magnetic field. Transitions from 
n = to final states of n > 2 have not been considered here 
because their matrix elements vanish in our approximation 
k±r± = 0. 

The procedures used here are much inferior to those for 
zero-field cross-sections as in the work of Reilman & Man- 
son (1979), in particular, the unperturbed final-state wave 
function and the state of the initial bound electron. The 
latter has been treated in the adiabatic approximation as 
a single-particle state with good quantum numbers n = 
and I. Use of this approximation in the case of atomic hy- 
drogen has been extensively studied by Potekhin, Pavlov & 
Ventura (1997). The true ground state contains small com- 
ponents with n > which lead to a non-zero value of the 
cross-section ctq which increases in significance as the thresh- 
old for the n = 1 final state is approached. The same is 
true for any single-particle state in a multi-electron ion. But 
in multi-electron ions, there are further problems. Whilst a 
shell model is clearly the only useful basis for description of 
the atom, as in the zero-field case, there must be a residual 
electron-electron interaction whose matrix elements, though 
satisfying the 5{n — I) = selection rule overall, can intro- 
duce an n = 1 and I + 1 component into the n = and / 
ground state and partially invalidate the adiabatic assump- 
tion. Relative to unity, the amplitude of such a component 
would be of order A/Hlob, where A is the matrix element, 
and is almost certainly very small. But the component has 
a photo-electric transition matrix element to the n = final 
electron state that is not inhibited by the sin x factor. Thus 
there must be a small non-zero value of ctq and a correction 
to Ctq. Clearly the cross-sections obtained here are subject 
to some uncertainty for the many reasons that have been 
discussed. However, we do not need great accuracy in the 
present work because the photo-electric transition rates re- 
quired arc also dependent on the high-frequency exponential 
tail of the blackbody spectrum and on the Lorcntz transfor- 
mation to the ion rest frame. But it is important that cross- 
sections with the correct high-field Landau level structure 
and dependence on photon polarization be used. 
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